Conformal metrics and true "gradient flows" for curves

We wish to endow the manifold M of smooth curves in Rn with a Riemannian metric that allows us to treat continuous morphs (homotopies) between two curves c0 and c1 as trajectories with computable lengths which are independent of the parameterization or representation of the two curves (and the curves making up the morph between them). We may then define the distance between the two curves using the trajectory of minimal length (geodesic) between them, assuming such a minimizing trajectory exists. At first we attempt to utilize the metric structure implied rather unanimously by the past twenty years or so of shape optimization literature in computer vision. This metric arises as the unique metric which validates the common references to a wide variety of contour evolution models in the literature as "gradient flows" to various formulated energy functionals. Surprisingly, this implied metric yields a pathological and useless notion of distance between curves. In this paper, we show how this metric can be minimally modified using conformal factors that depend upon a curve's total arclength. A nice property of these new conformal metrics is that all active contour models that have been called "gradient flows" in the past will constitute true gradient flows with respect to these new metrics under specific time reparameterizations

A Surface Evolution Approach to Probabilistic Space Carving

We present a 3D photography method that generates a texture-mapped three-dimensional model of a scene computed from multi-view calibrated two-dimensional photographs. Our approach first performs probabilistic space carving, which results in a 3D grid of voxel probabilities that describe the likelihood of a voxel existing in the model. We then employ a three-dimensional geodesic active surface to extract a most likely conformally-weighted minimal surface from the voxel probabilities. This surface is then polygonalized and texture-mapped, yielding in a 3D model easily rendered with standard graphics hardware

Deep learning-based methods for prostate segmentation in magnetic resonance imaging

Magnetic Resonance Imaging-based prostate segmentation is an essential task for adaptive radiotherapy and for radiomics studies whose purpose is to identify associations between imaging features and patient outcomes. Because manual delineation is a time-consuming task, we present three deep-learning (DL) approaches, namely UNet, efficient neural network (ENet), and efficient residual factorized convNet (ERFNet), whose aim is to tackle the fully-automated, real-time, and 3D delineation process of the prostate gland on T2-weighted MRI. While UNet is used in many biomedical image delineation applications, ENet and ERFNet are mainly applied in self-driving cars to compensate for limited hardware availability while still achieving accurate segmentation. We apply these models to a limited set of 85 manual prostate segmentations using the k-fold validation strategy and the Tversky loss function and we compare their results. We find that ENet and UNet are more accurate than ERFNet, with ENet much faster than UNet. Specifically, ENet obtains a dice similarity coefficient of 90.89% and a segmentation time of about 6 s using central processing unit (CPU) hardware to simulate real clinical conditions where graphics processing unit (GPU) is not always available. In conclusion, ENet could be efficiently applied for prostate delineation even in small image training datasets with potential benefit for patient management personalization

Artificial Intelligence for Classifying the Relationship between Impacted Third Molar and Mandibular Canal on Panoramic Radiographs

The purpose of this investigation was to evaluate the diagnostic performance of two convolutional neural networks (CNNs), namely ResNet-152 and VGG-19, in analyzing, on panoramic images, the rapport that exists between the lower third molar (MM3) and the mandibular canal (MC), and to compare this performance with that of an inexperienced observer (a sixth year dental student). Utilizing the k-fold cross-validation technique, 142 MM3 images, cropped from 83 panoramic images, were split into 80% as training and validation data and 20% as test data. They were subsequently labeled by an experienced radiologist as the gold standard. In order to compare the diagnostic capabilities of CNN algorithms and the inexperienced observer, the diagnostic accuracy, sensitivity, specificity, and positive predictive value (PPV) were determined. ResNet-152 achieved a mean sensitivity, specificity, PPV, and accuracy, of 84.09%, 94.11%, 92.11%, and 88.86%, respectively. VGG-19 achieved 71.82%, 93.33%, 92.26%, and 85.28% regarding the aforementioned characteristics. The dental student's diagnostic performance was respectively 69.60%, 53.00%, 64.85%, and 62.53%. This work demonstrated the potential use of deep CNN architecture for the identification and evaluation of the contact between MM3 and MC in panoramic pictures. In addition, CNNs could be a useful tool to assist inexperienced observers in more accurately identifying contact relationships between MM3 and MC on panoramic images

Performance of radiomics features in the quantification of idiopathic pulmonary fibrosis from HRCT

Background: Our study assesses the diagnostic value of different features extracted from high resolution computed tomography (HRCT) images of patients with idiopathic pulmonary fibrosis. These features are investigated over a range of HRCT lung volume measurements (in Hounsfield Units) for which no prior study has yet been published. In particular, we provide a comparison of their diagnostic value at different Hounsfield Unit (HU) thresholds, including corresponding pulmonary functional tests. Methods: We consider thirty-two patients retrospectively for whom both HRCT examinations and spirometry tests were available. First, we analyse the HRCT histogram to extract quantitative lung fibrosis features. Next, we evaluate the relationship between pulmonary function and the HRCT features at selected HU thresholds, namely -200 HU, 0 HU, and +200 HU. We model the relationship using a Poisson approximation to identify the measure with the highest log-likelihood. Results: Our Poisson models reveal no difference at the -200 and 0 HU thresholds. However, inferential conclusions change at the +200 HU threshold. Among the HRCT features considered, the percentage of normally attenuated lung at -200 HU shows the most significant diagnostic utility. Conclusions: The percentage of normally attenuated lung can be used together with qualitative HRCT assessment and pulmonary function tests to enhance the idiopathic pulmonary fibrosis (IPF) diagnostic process

Improved CTA Coronary Segmentation with a Volume-Specific Intensity Threshold

State-of-the-art CTA imaging equipment has increased increased clinician's ability to make non-invasive diagnoses of coronary heart disease; however, an effective interpretation of the cardiac CTA becomes cumbersome due to large amount of imaged data. Intensity based background suppression is often used to enhance the coronary vasculature but setting a fixed threshold to discriminate coronaries from fatty muscles could be misleading due to non-homogeneous response of contrast medium in CTA volumes. In this work, we propose a volumespecific model of the contrast medium in the coronary segmentation process to improve the segmentation accuracy. The influence of the contrast medium in a CTA volume was modelled by approximating the intensity histogram of the descending aorta with Gaussian approximation. It should be noted that a significant variation in Gaussian mean for 12 CTA volumes validates the need of volume-wise exclusive intensity threshold for accurate coronary segmentation. Moreover, the effectiveness of the adaptive intensity threshold is illustrated with the help of qualitative and quantitative results

Fully Automated Electrophysiological Model Personalisation Framework from CT Imaging

International audienceThere has been a recent growing interest for cardiac computed tomography (CT) imaging in the electrophysiological community. This imaging modality indeed allows to locate and assess post-infarct scar heterogeneity, allowing to predict zones of abnormal electrical activity and even personalise EP models. To this end, most of the literature uses manually segmented CT images where one fundamental information is extracted, the myocardial wall thickness. In this paper, we evaluate the impact of using an automated deep learning (DL) methodology to segment the left ventricular wall and extract relevant scar information on the resulting personalised models. Using CT images from 8 patients that were not used during the DL training, we show that the automated segmentation is very similar to the manual one (median Dice score: 0.9). Thickness information obtained this way is also very close to the manual one (median difference: 0.7 mm). A wavefront propagation model personalisation framework based on this thickness information does not show relevant differences in its output (median difference in local activation time: 2 ms), proving its robustness. Bipolar electrograms, simulated through a novel approach, do not differ significantly between manual and automated segmentations (Pearson's r: 0.99)

3D global and regional patterns of human fetal subplate growth determined in utero

The waiting period of subplate evolution is a critical phase for the proper formation of neural connections in the brain. During this time, which corresponds to 15 to 24 postconceptual weeks (PCW) in the human fetus, thalamocortical and cortico-cortical afferents wait in and are in part guided by molecules embedded in the extracellular matrix of the subplate. Recent advances in fetal MRI techniques now allow us to study the developing brain anatomy in 3D from in utero imaging. We describe a reliable segmentation protocol to delineate the boundaries of the subplate from T2-W MRI. The reliability of the protocol was evaluated in terms of intra-rater reproducibility on a subset of the subjects. We also present the first 3D quantitative analyses of temporal changes in subplate volume, thickness, and contrast from 18 to 24 PCW. Our analysis shows that firstly, global subplate volume increases in proportion with the supratentorial volume; the subplate remained approximately one-third of supratentorial volume. Secondly, we found both global and regional growth in subplate thickness and a linear increase in the median and maximum subplate thickness through the waiting period. Furthermore, we found that posterior regions—specifically the occipital pole, ventral occipito-temporal region, and planum temporale—of the developing brain underwent the most statistically significant increases in subplate thickness. During this period, the thickest region was the developing somatosensory/motor cortex. The subplate growth patterns reported here may be used as a baseline for comparison to abnormal fetal brain development